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This tentative correlation needs confirmation and explana- 
tion. If the jet regions of He 1 1.0830yum emission are closely 
associated with those producing X-rays, high resolution infared 
spectroscopy can be used to observationally study the velocity 
fields in the hot plasma regions of the jets. This would pro- 
vide the necessary evidence to test and further develop the- 
oretical models of intermediately fast (> 500 - 1500km s _1 ) 
interstellar shock waves. The HH 154 jet flow from the em- 
bedded protostellar binary L 1551 IRS 5 provides a case study, 
since adequate IR and X-ray spectroscopic data are in exis- 
tence. The thermal X-ray spectrum is fed into a photoionization 
code to compute, in particular, the line emission of He I and H I 
and to account for the observed unusual line intensity ratios. 
The advanced model is capable of accounting for most observ- 
ables, but shows also major weaknesses. It seems not unlikely 



, that these could, in principle, be overcome by a time depen- 
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dent hydrodynamical calculation with self-consistent cooling. 
However, such sophisticated model development is decisively 
beyond the scope of the present work. Continued X-ray obser- 
vations, coordinated with simultaneous high resolution infrared 
spectroscopy, are highly desirable. 

Key words. ISM: Herbig-Haro objects - jets and outflows - 
individual objects: L 1551 IRS 5 - Stars: formation - pre-main 
sequence 
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Abstract. The high energies of protostellar jets, implied by recent observations of X-rays from such flows, came very much as 
a surprise. Inferred shock velocities are considerably higher than what was previously known, hence putting even larger energy 
demands on the driving sources of the jets. The statistics of X-ray emitting jets are still poor, yet a few cases exist which seem 
to imply a correlation between the presence of He 1 1 .0830 /mi emission and X-ray radiation in a given source. 



1. Introduction 

The radiation from HH-objects and collimated jet flows 
is a manifestation of interstellar shock waves, the typical 
velocities of w h ich we r e believed t o be V., < s 500 km s 1 (e.g., 
Hartigan et all Il987l iRaga et all Il996t iReipurth & BalM 
200*lF Unexpectedly, a number of such objects from 
embedded prot ostars have r e cently b een detected at X- 
ray energies JPraydo et all 1200 it fravataetall |2002[ 
Preibischl 120031 IPravdoetall kool Tsuiimoto et al., 
20041 Idiideletall 120051: Idrossoetall EoOfil with refer- 
ence to HH 1/2, HH 1 54, Ser SMM 1 , HH 80/8 1 , TKH 8, HH 158, 
HH210, respectively). 

The X-ray spectra have typically a thermal spectral dis- 
tribution, characteristic of temperatures of several million de- 
grees. This was very surprising, since implied shock velocities 
(» 300 km s~ 1 ) would be very much higher than what had pre- 
viously been inferred from spectroscopic observations at longer 
wavelengths. This indicated that highly significant jet material 
had escaped detection. As a consequence, there is a potential 
risk that energy requirements on the jet source, being propor- 
tional to the square of the jet velocity, had been seriously un- 
derestimated. Proposed energy extraction mechanisms from the 
protostar-disk system might need to be revised in the light of 
these new findings. 

As in active stars, the continuous X-ray emissio n may 
be ac companied by Hel 1.0830/mi line radiation (e.g., Zirin, 
Il982l) . If that is the case, this line would be a valuable tool to 
probe the velocity fields of the hottest plasma regions in the 
jet. Also, in the infrared, extinction problems would be con- 
siderably reduced compared to shorter wavelengths and, for 
consistency checks, recombination lines from hydrogen, i.e. 
P/J 1 .28 18 //m and Py 1 .0938 /mi, could also be observed simul- 
taneously in the photometric J-band. 

In Tabled observations of Hel 1.0830//m and X-rays to- 
ward various flows have been compiled. Evidently, many ob- 
jects lack the Hel 1.0830/mi line, but reveal a rich H2 spec- 



trum, which is consistent with non-dissociative shocks of low 
velocity, of the order of lOkms 1 . This is true also for HH 1G 
and the SerSMMl jet, toward which both Hel 1.0830/mi and 
X-ray observations have been performed, but neither of which 
was detected. HH 2H was detected in X-rays but has apparently 
not been observed spectroscopically in the 1 /mi region. On 
the other hand, the jet flows HH 158 (DGTau A) and HH 154 
(L1551 IRS 5) were clearly detected in both Hell.0830/mi 
and at X-ray wavelengths, lending support to the above propo- 
sition. 

It is thus the primary objective of this Research Note to 
stimulate further X-ray observations for the investigation of the 
high energy regions in protostellar jets, alongside with coordi- 
nated high resolution near infrared spectroscopy to determine 
Doppler velocities and line shapes. These data can then be com- 
pared with predictions from theoretical shock models. 

To date, the best available data are those of the HH 154 
flow and we shall first examine the accumulated evidence for 
this object (Sect. 2). Some numerical model calculations will 
be presented and discussed in Sect. 3. In that section, we will 
also conclude with some suggestions for future work. 

2. L1551 jet (HH 154): the data 

2.1. Infrared spectroscopy 

The exami ned data are bas ed on the J- and H-band spec- 
troscopy bv lltoh et alJ (J2000). Forbidden lines of ionized iron, 
[Fe II], dominate the observed s pectru m and this has been anal- 
ysed in depth by iLiseau et all d2005l) . The spectral resolution 
near A = 1 .0 /mi is quite modest, AA = 3.6x1 0~ 3 /mi (R = 280, 
AV = 10 3 kms _1 ), which also corresponds to the accuracy 
with which wave lengths can be read off Fig. 5 in the paper by 
lltohet alJ J2000l) . 

In Table|2j the line fluxes are given and in Fig.^ the jet- 
spectrum between 1.05 /mi and 1.35 //m is displayed. Include d 
is the best fit for the [Fe II] transitions of Lisea u et alJ 
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Table 1. Protostellar jets with X-ray data and/or NIR spectra 



Jet flow 


He 1 1.0830 /im 


X-rays 


Reference 


Comment 


HH 158 






1 2 


TFp TT1 
t rc 11 J 


HH 154 


v 

A 


y 


1 d < 


TFp TT1 
[re ii j 


HH 210 




y 






HH 1G 


A 


y 

A 


7 8 
/ . o 


m i [rc ii j 


HH 2H 




v 

A 



o 




TKH 8 




A 


Q 




HH 43 






10 


112 


HH 24A 


_ y 




1 1 


Hi + TFe TT1 
in i [rc ii j 


[111 Z.J 


— X 




1 1 

1 1 


H, 


HH 26 


_ v 
A 




1 1 




HH 72 


— X 




1 1 


H, 


VelIRS17 


-X 




12 


H; 


HH320 


-X 




11 


H; 


HH 321 


-x 




11 


H 2 


HH 80 




X 


13 




HH81 




X 


13 




SerSMMl 


-X 


-X 


14, 15 


H 2 


HH 99-BO 


X 




16 


H 2 + [Fell] 



Notes to the table: x = detected, -X means observed but not detected. 
H 2 indicates the presence of ro-vibrationally excited H 2 lines; [Fell] 
the presence of [Fell] lines; H 2 + [Fell] the presence of both species 
along the same line of sight/insi de spectrograph slit. When possible, 



HH-nomenclatur e follows that ofjRe jurtbl il999j)__ 
Reference s : (1) iTakami et all <200?k (2 1 lOudel et all 1200 




shown by the smooth red line, whereas the obs ervations are 
shown as black histograms (see also Tabled iLiseau et alJ 
J2005I) tentatively identified the feature near 1.08 /mi as He I 
3 P" - 3 Si 1.083034/mi, and in the following, possible other 
identifications, because of wavelength overlaps within the er- 
rors, are briefly discussed. 



2.1.1. [Fell]1.0890/im 

A little red "bump" is visible under the observed line near 
1 .08 /im. This bump is due to [Fe II] a 4 H - b 2 G at 1 .08899 jum, 
with upper level energy E up /k > 44 000 K. The matching of 
the observed line with this transition would require, therefore, 
a much higher level of excitation than what has been inferred 
from the other lines in the spectrum, and is therefore inconsis- 
tent with the observations. Consequently, [Fell] 1.0890/mi is 
disregarded as a viable candidate for the 1.08 //m feature. 

2.1.2. [SI] 1.0821 /im 

Around 1 .08 //m, atomic/ionic ground state transitions exist 
only for neutral sulphur. Relatively close in wavelength to the 
observed feature is [SI] *D 2 - 3 P 2 at 1.082117/mi. The other 
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4x10-12 



2xl0-' a 



c 
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3 



L1551 IRS 5 Jet-N 



CD P- 
O CO 
CO 10 




rms noise level 




1.2 

Wavelength, A(/zm) 



1.3 



Fig. 1. The J-band spectrum of lltoh et all J2000l) containing 
the Paschen lines, with the estimated level of the rms noise 
indicated. The observations are displayed as the black his- 
togram, whereas the smooth red line refers to the fit of 
the [ Fe II] spectrum (for Ay = 5.7 mag, ILiseau et all 
120051) . The wavelength positions of the recombination lines 
Hell.0830/mi, Py 1.0938 /mi and P/3 1.28 18 /mi are indicated 
above and below the spectrum. This applies also to the wave- 
lengths of the transitions to the inverted ground state of [S I], 
i.e. l D 2 - 3 P 2 1.08212/mi, l D 2 - 3 Pj 1.13058/mi and 'D 2 - 3 P 
at 1.15376/im, respectively. 



two transitions to the inverted ground state of sulphur are also 
in the J-band, viz. [ D 2 - 3 Pi at 1.130585//m and 'D 2 - 3 P 
at 1.1537564yum (Fig.Q. None of these latter lines are de- 
tected. This is hardly surprising, as sulphur is fully ionized 
in the jet (j^si = 10.36 eV), which is further supported by the 
non-detection of the [S I] 25.25 yum fine structure 3 Pi - 3 P 2 line 
dWhite et all l200oT) . The observed line near 1.08 /mi is there- 
fore not due to neutral sulphur. 



2.1.3. Pyl.0938/im 

Strong B ai mer line emission ha s frequently been observed 
from the jet llFridlund et alll2005[ and references therein) and, 
in view of the high extinction, Paschen lines could be expected 
to be detectab le in the J-band. For insta nce, using the extinc- 
tion curve of iRieke & LeboTskl lll985h for the value of Av 
= 5.7 mag (derived from the [Fell] spectrum), the attenuation 
('reddening') of the Ha flux corresponds to a factor of 70. In 
contrast, at the wavelength of P/3, this extinction amounts to a 
diminishing factor of only 4. However, based on Ho- observa- 
tions of the jet, the 2.5<r feature at 1.2785/mi would appear to 
be too strong (if real) by a factor of about 3 to be consistent with 
P/3 1 .28 18 //m. As is evident from Fig.^ this feature is instead 
well fit by an [Fe II] line. 

For a wide range in density, intrinsic Case B line ratios are 
typically P/3/Py ~ 2 jHummer & Storevl Il987l) . Since the ob- 
served flux ratio F(1.2785)/F( 1.0866) <K 1, and any extinc- 
tion correction could only decrease this value even further, the 
observed feature near 1.08 /mi is not due to P-y. 
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Table 2. Wavelength and flux estimates in t he J- and H-band 
spectra of HH 154 (based on observations bv lltoh etal 



/l obs 

(/an) 



10 15 X Fy 

line 

(erg cm -2 s _1 ) 



Identification & A.^ (jim) 
(Fell: Ouinet et al.. 19961 



1.0866 


26.0 ± 4.6 


He 1 1.083034 3 


P°2p 


- 3 S2s 


1.2566 


70.2 ± 4.6 


[Fell] 1.256680 


a 4 D 


-a 6 D 


1.2785 


13.2 + 5.2 


[Fell] 1.278776 


a 4 D 


-a s D 


1.2942 


19.2+ 1.6 


[Fell] 1.294268 


a 4 D 


-a 6 D 


1.3205 


25.5 + 2.1 


[Fell] 1.320554 


a 4 D 


-a 6 D 


1.3272 


19.9 + 4.0 


[Fell] 1.327776 


a 4 D 


-a 6 D 


1.5334 


20.7 +1.1 


rr ' TTi 1 rio i 

[Fell] 1.533471 


a 4 D 


- a 4 F 


1.5995 


13.0 + 1.7 


TFeTTl 1 S99473 


a 4 D 


- a 4 F 


1.6435 


112.3 + 1.9 


[Fell] 1.643550 


a 4 D 


-a 4 F 


1.6678 


8.6 + 2.0 


[Fell] 1.663766 


a 4 D 


-a 4 F 


1.6769 


13.7 + 1.7 


[Fell] 1.676876 


a 4 D 


-a 4 F 


1.7045 


5.0 + 2.0 








1.7156 


6.0 + 2.0 


[Fell] 1.711129 


a 4 D 


-a 4 F 


1.7490 


6.0 + 2.0 


[Fell] 1.744934 


a 4 D 


-a 4 F 


1.8045 


13.0 + 4.0 


[Fell] 1.800016 


a 4 D 


-a 4 F 


1.8134 


22.0 + 6.6 


[Fell] 1.809395 


a 4 D 


-a 4 F 


1.0938 


1.4 ±4.6 


HIPy 






1.2818 


4.1 +5.2 


HIP/3 







Notes to the Table: 

Estimated accuracy is AA = ±0.004 yum. 
* Flux calibration from Y. Itoh (private communication). 



2.1.4. He 1 1.0830 

For reasons which will become apparent below, we find 
that He 1 1.0830/im is the best identification for the observed 
1.08 yum feature. 

Successful models of the infrared [Fell] spectrum of the 
jet implied Ay = 5.7 mag, where an extinction law of the form 
A,,/A v = 3.0134 - 5.62U + 4.3847 A 2 - 1. 5839A 3 + 0.2181A 4 . 
0.55 < A < 2.2 //m had been adopted ( Ri eke & Lebofskvl 
Il985h , When the extinction is spectroscopically determined, 
the uncertainty in Ay is often dominated not by the measure- 
ment errors, but by the uncertainty in the atomic parameters for 
the transitions from common upper levels (Tableland Fig. 1 of 
Liseau et al., 2005). For instance, the widely exploited line ra- 
tio F( 1 .64 ,um) / F( 1 .25 ;um) is not be tter determined than to the 
range 0.74-0.96 JOuinet et allll996l) . where the higher value is 
based on more recen t calcu lations, and which was the preferred 
one bv lLiseau et alJ J2005h . For the observed ratio, extinction 
estimates would then fall into the range Ay = 5.7 - 8.6 mag 
(subject to the assumed extinction curve). 

The extinction corrected luminosity of the He 1 1.0830/im 
line is Li.o830/jm = 6 x 10 28 erg s _1 (D = 140pc and assum- 
ing isotropic emission). This corresponds to 20% of the X-ray 
luminosity from this region in the jet (see next section). 



2.2. X-rays from the 0.3keV- 7.9 keV band 

At the distan ce of 140 pc , the X-ray luminos i ty is ( 2 - 5) x 
10 29 erg s _1 dFavata et all 120021 bonito et all |2004 . As can 
be seen in Fig.|3 these XMM-Newton data are consistent 
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Fig. 2. The observed X-ray spectrum of the L 1551 IRS 5 
jet is shown by the big dots, omitting error bars for clar- 
ity. The data are limited to roughly 0.3 keV - 3.3 keV. 
The spectrum appears thermal and is consistent with tem- 
peratures of about 4 - 4.5 MK (including obs ervational er- 
rors, t he best fit MEKAL model yields 4+2.5 MK. lFavata et all 
120021) . For the model computations, the spectral shape of the 
blackbody, shown by the continuous curves, has been used as 
input spectrum (see the text). The frequencies corresponding to 
the ionization edges of He I and He II are shown by the vertical 
dashed lines. 

with a thermal sou rce at temperatures of roughly 4 x 10 6 K. 
iFavata et all d2002h give the emission measure of the hot gas as 
EM = J(x e nu) 2 dV = 1.1 x 10 52 cirT 3 , yielding an estimate of 
the source size, the linear scale of which corresponds to a few 
xlO AU for densities of the order of 10 4 cm 4 . 

A large body of observationa l data exi s ts for t he X-ray spot 
within the (northern) jet jFridlund et all l2005t iLiseau et all 
120051) . The value of the X-rav attenuation. Av = (7.3+2.1)mag, 
is consistent with the visual extinction derived from the in- 
frared [Fell] spectrum. Conforming with thi s steep extinc- 
tion gradient is also the observed Hc/H/J-ratio dFridlund et all 



120051 

The analysis of the X-ray data for protostellar jets (Table[Q 
has generally been in terms of shock waves and has, as such, 
frequently exploited analytical formulations for the tempera- 
ture and pre-shock density, with the shock velocity as a parame- 
ter (e.g.. lOstriker & McKeelll988HRaga et alll2002l) . Common 
results can be summarized as T > 10 6 K, V s > 300kms _1 and 
«o ~ a few x 10 2 crrT 3 (see the references in Sect. 1). 



3. Discussion and conclusion 

The overall observational evidence supports the view that the 
X-ray emission originates from fast shock s in the HH 154 jet, 
with velocities in excess of 500 km s 1 dFavata et all 120021 
iBonito et all 120041 Liseau et all Eool IFavata et all l2006l)~ 
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Fig. 3. Energy level diagramme of He I. For clarity, levels for 
quantum numbers > 4 have been omitted. The transitions, 
which are discussed in the text, i.e. 1.0830pm and 0.5876pm 
of the triplet states and 2.06 pm of the singlets, are marked and 
the wavelength of He I ionization is shown next to the slanted 
dashed line (see also Fig.|2}. 



reasonable scenario features a hypersonic jet, ramming into 
molecular cloud material which leads to the complete destruc- 
tion of molecules and grains and the ionization of the gas, 
which is heated to high temperatures. This hot plasma is ra- 
diatively cooled by X-rays, which are capable of ionizing the 
neutral gas ahead of the shock, altering its state of excitation 
and potentially also its chemistry. In order to arrive at quantita- 
tive estimates of the associated helium and hydrogen recombi- 
nation line emission, in particular, we have run photoionization 
models for this scenario. 

The He 1 1.0830 pm line originates from levels which are 
more than 20 eV above the ground (Fig. QJ. Evidently, the ob- 
served X-rays are sufficientally energetic to photoionize helium 
(see Fig.|2ji resulting in He 1.0830 emission upon recombina- 
tion. But is it possible to explain the fact that the observed line 
intensity ratio He 1 1.0830pm/Py 1.0938 /an » 1? Such in- 
verte d ratios are gen e rally n ot revealed by models of jet shocks 
(e.g. lHartigan et all Il987l) . but have been observed to de- 
velop in supernova ejecta fe.g. lMeikle et all|l993HPozzo et all 
12004 . 

For the model calculatio ns we used the prog ram Cloudy 
(version 96 last described by iFerland et all Il998l) . which de- 
termines the temperature structure for a gas of a certain ge- 
ometry and density, including the transfer of radiation. For the 
calculations described here, the shape of the X-ray input spec- 
trum is that of a blackbody (Fig.|2j, scaled by the observed 
X-ray luminosity, i.e. the inner radius is r- m 
The outer radius of the spherical cloud was initially set by the 
observed X-ray emission measure and the assumed constant 



density, i.e. r ollt ~ (3 EM/4n) l/i n(H) 2/3 . Solar abundances 
were assumed for the chemical composition of the gas, hence 
«(He)/«(H) = 0.1. 

The resulting model for n(H) = 10 5 crrT 3 , which is in 
fair agreement with the observations, is depicted in Fig.@] 
In particular, the ratios He 1 1.0830pm/Py 1. 0938 pm and 
HeI1.0830pm/P/n.2818pm are 13 and 7, respectively, i.e. 
both much larger than unity as has been observed (Table[2). 
According to the model, with e -collision coefficients from 
iBrav et alJ d2000h . the excitation of the He 1 1.0830pm line is 
dominated by collisions, so that in this case, the '10830 re- 
combination line' is not really due to He recombination (com- 
pare the JCHeii and jfHei regions in the figure). For purely ther- 
mal broadening, the line would be marginally optically thick, 
i.e. r = 1.5 at the center of the line. The flux of the reddened 
2.06pm line (Fig.[5} is predicted not to exceed the 10% level 
of that in the 1.0830pm line. 

A hint to the nature of the He I excitation is potentially 
also provided by the He 1 0.5876pm line, originating from the 
state just above that of the Hel 1.0830pm line (see Fig.[3j. 
The line intensity ratio He 1 0.5876 pm/He 1 1.0830pm of the 
model is 0.3, i.e. the predicted (reddened) He 1 0.5876 pm flux 
is 4 x 10~ 16 erg cm' 2 s" 1 . This is i n agreement with the ob- 
servation by ICohen & Fuller! dl985l) . which resulted in an up- 
per limit, i.e. F < 6.5 x 10~ 16 erg cirT 2 s _1 (the units in that 
paper seem erroneous ). Given the information provided by 
ICohen & Fuller! dl985l) . it seems not likely, however, that this 
measurement actually refers to the location of the X-ray source. 
In addition, the X-ray emission is known to be variable on short 
time scales dFavata et alll2006l) . 

The observed Hel 1.0830pm luminosity (Sect. 2.1.4) is a 
sizable fraction of the X-ray luminosity and it comes hardly as a 
surprise that the model falls short with regard to Li.o830^m- This 
may not be attributable solely to source variability. In principle, 
an ad hoc increase of the emitting volume could artificially cure 
this luminosity problem: for instance, shock models would lead 
to typically much larger cooling lengths (10 15 - 10 16 cm) than 
what is indicated in Fig.|4](< 10 14 cm). The required increase 
by an order of magnitude seems consistent with th e models 
of X-ray emitting bow shocks of iRaga et all d2002l) . Also, as 
suggested by the referee, a jet with a time variability in the 
ejection velocity could increase the emitting volume (hence the 
total luminosity). This volume increase is produced when two 
successive working surfaces collide with each other. 

However, such an ad hoc procedure would lack any physi- 
cal justification. Therefore, rather than trying 'to fix' this poor 
model, one would like to see a proper time-dependent hydro- 
dynamical calculation for the high velocity jets, accounting for 
the shocks and the cooling in a self-consistent manner. This is 
also warranted by the various time scales of the steady state 
photoionization model (e.g., for thermal equilibrium and hy- 
drogen recombination), which are all very much longer than 
the dynamical time scale of less than one year. 

On the observational side, one would like to see concerted 
efforts of X-ray observations and simultaneous intermediate to 
high resolution (R > a few x 10 3 ) infrared spectroscopy. 
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Fig. 4. Helium ionization structure of gas with constant density 
«(H) = 10 5 cirT 3 , being illuminated by a 4MK point source. 
The run of the temperature T e (in units of 10 6 K) is shown in 
red, that of the electron fraction x e in green. The fractional 
H I density is shown by the black broken line, whereas those 
of He I, He II and He III in blue, cyan and magenta, respec- 
tively. Tests with the previous version Cloudy 90 gave essen- 
tially identical results, except that the temperature structure is 
smoother. The hump due to 'ionization jumps' for a number of 
species is absent, indicating that this feature does not have any 
major effect on the results of interest. 
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